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ABSTRACT 

Aims. To assess the nature of three WNh stars in the SMC, to constrain stellar evolution beyond the main sequence at low metallicity, 
and to investigate the metallicity dependence of the clumping properties of massive stars. 

Methods. We compute atmosphere models to derive the stellar and wind properties of the three WNh targets. A FUV/UV/optical/near- 
infrared analysis gives access to temperatures, luminosities, mass loss rates, terminal velocities and stellar abundances. 
Results. All stars still have a large hydrogen mass fraction in their atmosphere, and show clear signs of CNO processing in their 
surface abundances. One of the targets can be accounted for by normal stellar evolution. It is a star with initial mass around 40-50 
M G in, or close to, the core He burning phase. The other two objects must follow a peculiar evolution, governed by fast rotation. 
In particular, one object is likely evolving homogeneously due to its position blue-ward of the main sequence and its high H mass 
fraction. The clumping factor of one star is found to be 0.15+0.05. This is comparable to values found for Galactic Wolf-Rayet stars, 
indicating that within the uncertainties, the clumping factor does not seem to depend on metallicity. 
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1. Introduction 

Massive stars have gained considerable interest in the last years 
in various fields of astrophysics. In addition to being well rec- 
ognized actors in Galactic astronomy (HII regions, bubbles, 
ISM dynamics and winds, triggered star formation, chemical 
enrichment) their role in extragalactic physics has considerably 
grown. Massive stars are kn own to govern starb urst events lo- 
cally and at high re dshift (e.g.lSteidel et alj [l996). being respon- 
sible for outflows dVeilleux et alj 120051) and inte nse episodes 
of sta r formation in luminous infrared galaxies dGenzel et alj 
1 1998b. The association of some GRBs with c ore collapse super- 
nova dGalama et al.ll 1 998t [Hiorth et al.ll2003l) make fast rotatin g 
massive stars the likely progenitors of GRBs dWooslevlll993l) . 
However, the conditions under which a massive star will actu- 
ally end up in such an event are still poorly known. Fast rotation 
is of course a key parameter, but the way the star loses angu- 
lar momentum, and thus its mass loss history, is equally im- 
portant. With mass loss decreasing at lower metallicity, GRBs 
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are thought to occur in low Z environments. Observations of 
GRB hosts tend to confirm this expectation dProchaska et alj 
2004). However, the existence of an upper metallicity limit 
for the formation of GRBs is not established. Massive stars 
are also thought to be responsible for the re-ionization of the 
Universe at very high redshift (e.g. Schaerer at al. 2002 and 
references therein). Indeed, the first (population III) stars may 
have been very massive objects with very hi gh effective temper- 
atures due to their metal free composition (Bro mm et al.l ll999: 
iNakamura & Umemural 19991) . 

In spite of this growing importance of massive stars, their 
evolution is not fully understood even in the local Universe. It is 
not clear at present which phases a star of a given initial mass 
will go through during its life. For example, the initial condi- 
tions for a star to enter the red supergiant or luminous blue vari- 
able stages are rather blurry. Similarly, the evolutionary connec- 
tions between various types of Wolf-Rayet stars is poorly con- 
strained. These difficulties result from two facts. The first one 
is that massive star evolution is governed by mass loss, and our 
knowledge of mass loss rates through the upper part of the HR 
diagram is still patchy. The second problem is due to rotation: 
stars with similar initial masses but different rotational veloci- 
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ties will follow different paths in the HR diagram, and th us will 
evolve through different states dMevnet & Maeder 2000). These 
two ingredients (mass loss and rotation) are the keys to mas- 
sive star evolution. An additional complication occurs because 
both of the above factor s depend on m etallicity: stars lose less 
mass at low meta llicity (Puis e t al.| 20Q 0j), and at the same time 
they rotate faster (Marta van et al.ll2007l) . Consequently, it is not 
only crucial to know the mass loss and rotation of all types of 
massive stars, but it is also important to constrain their metal- 
licity depen dence. Recent results i ndicate a M ~ Z° 83 behavior 
for O stars dMokiem et al.ll2007bl) . Theoretical predictions and 
observation al constraints also exist for various types of Wolf- 
Rayet stars (IVink & de Koterll2005t iGrafener & Hamannll2008t 
ICrowther et alj 2002al) . A rather clear picture of the relative 
strength of stellar winds at different metallicities is thus emerg- 
ing. But the question of their absolute strength is not settled. The 
reason is clumping. 

There is now no doubt that massive-star winds are not ho- 
mogeneous. Evidence comes from both indirect and direct trac- 
ers. The f ormer compri se electron scattering wings of emis- 
sion lines ( Hiilierlll99ll). too-weak observed UV line profiles 
(ICrowther et alj I2002U iHillier et al.1 l200l iMassa et all l2003t 
IBouret et al.ll2003l 120051: IFullerton et al.ll2006l) . shapes of X-ray 
lines (Oskinovaet al. 2004), submm excess in continuum emis- 
sion (Blo mme et al.ll2002l), linear-polarization continuum varia- 
tions dSt.-Louis et al.lfl993l) . Direct indications are the observa- 
tions of global moving substructures on emission-l ine profiles 
dEversberg et al.lll998tlLepine & Moffa1l999ll2008h . Clumping 
affects mass-loss determination of massive stars because it 
changes the density, which in turn modifies the shape of wind 
lines. In practice, clumping is usually quantified by means of a 
filling factor /. For a given /, recombination lines will be fit- 
ted, to first order, with a mass loss rate approximately 1 / JJ 
times lower compared to a homogeneous atmosphere model. 
Knowledge of the clumping properties of massive-star winds 
is thus a key to the quantitative understanding of massive star 
evolution. 

Studies of Galactic stars indica te f f a ctors ranging from 
0.01 to 0.2 dHamann & Koesterkd 119981: Iffillier et : alj l200lt 
IBouret et al.ll2005t IFullerton et al. 2006]). Such facto rs may not 
be co nstant throughout the wind of a massive star (Pui s et al.l 
120061) . In addition, if inh omogeneities are due to line driving in- 
stabilities dRunacres & Owocklll2002l) . a metallicity dependence 
can be expected. It is however not c onstrained at present. A 
first step in this direction was made by iMarchenko et ail d2007) 
(hereafter paper I). In this study, the dynamical properties of 
clumps directly observed in emission lines of three SMC WR 
stars were found to be very similar to those of Galactic objects. 

In the present study, we analyze the three stars observed in 
paper I by means of atmosphere models. We constrain the stellar 
and wind properties, which help us understand their evolution- 
ary status. Being located in the SMC, these stars give us a clearer 
picture of stellar evolution beyond the main sequence in a low 
metallicity environment. We find that evolution at very high ro- 
tational rate is most likely at work, as suspected for the formation 
of GRBs. In addition, one can pinpoint the value of the clump- 
ing factor / for one star, showing no systematic difference with 
values found for Galactic stars. 

In the following, we briefly review the observational material 
in Sect. |2]before explaining the modelling technique in Sect. [3] 
We then present the results for each star in Sect. [4] and discuss 
them in Sect. [5] Sect.|6]summarizes our conclusions. 



2. Spectroscopic data 

The three stars analyzed here (SMC-WR1, SMC-WR2, SMC- 
WR4) were observed with the UV- Visual Echelle Spectrograph 
(UVES) on VLT/UT2 at ESO. Observations were conducted on 
August 27-28"' 2006. The regions between 3930-6030 was tar- 
geted. A slit width of 1.5" was used and observations were con- 
ducted under an airmass varying between 1.5 and 2.7, and a 
seeing fluctuating between 0.8 and 2". The data were reduced 
using the ESO UVES Pipeline. The reduction process includes 
bias and inter-order background subtraction, for both the objects 
frames and the flat-fields calibration images, optimal extraction 
of the objects (that includes removing skylines and cosmics re- 
jection), division by a flat-field extracted with the same weighed 
profile as the object, wavelength calibration and merging of all 
overlapping orders after a rebinning to a constant wavelength. 
No spectrophotometric standard was observed, so our spectra are 
not flux calibrated. The final spectra have a resolving power of 
~ 20000 and a signal to noise ratio between 120 and 160. More 
d etails about the observa tions and data reduction can be found 
in IMarchenko e t al. (2007). The spectra we use in the following 
are the combination of all spectra obtained during the observing 
run. 

We have also used the weighted mean sp ectra covering 
the H a region presented by Foell mi et al.l d2003l) . They have a 
low resolution (R~800) and result from observations with var- 
ious telescopes an d spectrographs (see details in Table 2 of 
iFoellmi et al.l l2003). These spectra are not flux calibrated. 

In addition to the optical spectra, we used archival IUE data 
for stars SMC-WR2 and SMC-WR4 (spectra SWP7026 and a 
weighted mean - according to exposure time - of SWP06195 
and SWP10727 respectively). They were taken in the low res- 
olution mode (R ~ 250) of the instrument. Two high resolution 
spectra exist for SMC-WR4, but their S/N is low and they clearly 
suffer from reduction artifacts (e.g. oscillating continuum). For 
those two stars we also retrieved FUSE archival data (spectra 
G036101000 and G03900301000). They have a resolving power 
of about 20000. 

Finally, ESO-SOFI spectra of the sample stars covering the 
J and H band (also K-band for SM C-WR4) were us ed. The data 
on SMC-WR4 were presented bv lCrowtherj d2000t) . For SMC- 
WR1 and SMC-WR2, we used unpublished data kindly provided 
by Paul Crowther. They were obtained in November 2003 (only 
SMC-WR1) and November 2004 (SMC-WR1 and SMC-WR2) 
with a 0.6" slit. Their resolution is ~ 500 and they are not flux 
calibrated. They were reduce d following the same procedure as 
described in lCrowfherl d2000l) . 



3. Modelling 

The stellar and wind properties of the stars analyzed here have 
been derived using non-LTE, line blanketed models includ- 
ing winds. The models have been computed with the code 
CMFGEN. A detailed description of the code is given in 
iHillier & M iller (1993). Here, we recall only a few points rel- 
evant for the understanding of the present study. 

Line-blanketing: a super-level approximation is used to al- 
low the inclusion of a large number of atomic levels from var- 
ious elements. In our models, H, He, C, N, O, Ne, P, Si, S, Fe, 
Ni have been included. In total, the populations of ~1 150 levels 
and ~3500 super-levels were calculated through the resolution of 
the statistical equilibrium equations. A constant microturbulent 
velocity (20km s~') was adopted for the calculation of the atmo- 
spheric structure. For the formal solution of the transfer equa- 
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tion leading to the emergent spectrum, a microturbulent velocity 
ranging from 10-50 km s _1 (depending on the star) to O.lXVoo 
from bottom to top of the atmospher e was chosen. Stellar abun - 
dances were chosen to be l/5' h solar dGrevesse & Sauva| [l998). 

Clumping: CMFGEN allows the inclusion of clumping with 
the following formalism. A volume filling factor / is assumed 
to vary monotonically with depth, starting from a value of 1 at 
the photosphere and reaching a maximum value f m at the outer 
boundary of the atmosphere. In practice, / is parametrized as a 
function of the local wind velocity v 



/ = foo + (1 - /«)« 



(1) 



where v c ; is a parameter fixing the velocity at which clumping 
starts to become significant. In practice, we chose v c i = 100 
km s _I so that clumping starts above the sonic point. L ower 
values have been reported for O stars (Bouret et al. 2005) but 
are difficu l t to ex plain theoretically (Runacres & Owockil l2002). 
iPuls et all d2006l) showed that in O stars, clumping varied non 
monotonically: after an increase (/ decreases) up to the region 
around the sonic point, clu mping decreases in the out er wind. 
Theoretical simulations by iRunacres & OwockH d2002l) tend to 
predict a similar trend. However, this behavior has not been con- 
strained for Wolf-Rayet stars. Hence, we adopt a conservative 
approach and keep on working with a monotonic clumping law, 
keeping in mind that it might not be fully realistic. 

Density structure: Until recently a major problem with the 
analysis of W-R stars was that the density structure had to 
be specified using an assumed ve locity law. However work by 
Graf ener & Hamannl d2005l 12008b have allowed velocity laws 
(and mass loss rates) to be derived for some W-R stars using 
a formulation of radiation driving applicable to optically thick 
winds. A significant limitation of the work is that the volume 
filling factor /, and its variation with radius, must be assumed 
since as we have seen clumping properties are still poorly con- 
strained. 

Unfortunately, CMFGEN is not yet able to self-consistently 
compute the velocity law. It still assumes some analytical 
expression above (approximately) the sonic point: the so-called 
ft velocity law is connected to the photospheric structure. 
Below the sonic point, a photospheric stru cture is either 
param etrized by a scale height formalism ( Hillier & Miller 
1998) or is taken from alternative atmosphere models, usually 
TLUSTY (L anz & Hubenvl 120031) . In the present analysis, 
we have improved this formalism. We now require that hy- 
drostatic equilibrium is satisfied below the sonic point. In 
practice, we proceeded as follows. We first started with a 
photospheric TLUSTY structure connected to a ft velocity 
law. With this structure, we started to converge the models. 
After a few iterations (usually 10-15), the predicted level 
populations were used to compute the radiative acceleration. 
This acceleration was subsequently injected in the hydrostatic 
equilibrium equation from which an altered density structure 
was derived at depth, below the sonic point. This ensures 
a higher degree of consistency between the hydrodynamic 
and atmospheric structures. Above V(sonic)/2, we adopted 
the usual beta velocity law formulation (see below). With 
this revised structure, the models were converged for another 
10-15 iterations before another modification of the photospheric 
structure was performed. We ran 3^1 of these structure iterations 
before converging the atmosphere models. In practice, the main 
inconsistency in the final density structure (lack of radiative 
acceleration) was found around the sonic point. This reflects 
the fact that we are specifying the wind velocity law rather than 



deriving it, that the critical point is close to the sonic point, 
and that the disagreement depends on the clumping law adopted. 

With these models, we have constrained the main stellar and 
wind parameters using the following procedure: 

• Temperature. As usual for Wolf-Rayet stars, we have de- 
rived both an effective temperature representing the radia- 
tion flux at a Rosseland optical depth of 2/3, and a tem- 
perature T» at t=20. Both values can be significantly dif- 
ferent in the case of a very dense wind for which t=2/3 is 
reached in a region where the atmosphere is no longer static. 
In practice, the temperature of the model was adjusted to cor- 
rectly reproduce the relative strength of Nra /L14635,4641, 
Nra /U45 10,45 15, Niv ,14058, Nv ^4604, 4620 and Nv 
/14944 (when present). Infrared He lines - Hei A\. 08 pm, 
He ii /1 1.1 7 pm- were also used and usually gave consistent 
results (see Sect.© . In addition for one star (SMC-WR4) we 
could use the classical He i/He n optical line ratio since the 
optical He i lines were strong enough. We estimate the uncer- 
tainty on the derived temperatures to be of the order of 3000 
K for SMC-WR2 and SMC-WR4, 7000 K for SMC-WR1 
(for which N in is not present and N iv is weak). 

• Luminosity and extinction. For stars SMC-WR2 and SMC- 
WR4, the luminosity and local extinction was determined 
from the fit of the SED. The IUE (flux calibrated) spectra 
have been used in conjunction with optical photometry (our 
optical spectra are not flux calibrated) to adjust the stellar 
luminosity and extinction. Table Q] summarizes the photo- 
metric data for each star. FUSE data were not used since 
1 ) the extinction laws in the far-UV are not well constrained 
and 2) the IUE and FUSE spectra did not show the same 
flux level in the region around 1200 A (especially for SMC- 
WR2), pointing to an uncertainty in the flux calibration (ei- 
ther for FUSE or IUE data). Given point 1, we relied prefer- 
entially on IUE spectra. Infrared data were not used since the 
extinction law in the near-IR is not constr ained for the SMC. 
The Galactic (SMC) extinctions laws of lSeatonl d 19791) and 
iHowarthl d!983l) dBouchet et al.ll 1985b were used. The SMC 
extinction is stronger than the Galactic one in the UV (they 
are similar in the optical), so that using both IUE and opti- 
cal data helped constrain the local, SMC extinction and the 
foreground, Galactic E(B-V). A distance modulus of 19.0 
(corresponding to a dista nce of 63.2 kpc) w as assumed. This 
value was reported by di Benedetto] d 19971) and is interme- 
diate between low values (18.80) and high (19.20) values 



reported in the literature dCold 


19981 ICaouto et al.1 


1999 


Groenewegen 2000; Harries et al. 


2003; Hilditch etalj 


2005 



iKeller & Woodl2006l) . For star SMC-WR1, in the absence of 
(F)UV data, we only used the absolute V magnitude. From 
the observed B-V and the theoreti cal value expected for hot 
massive stars ((B-V) =-0.28, see iMartins & Rlezl 120061) we 
derive E(B-V)=0.24 which, together with an assumed dis- 
tance modulus of 19.0, implies My=-4.57. We adjusted the 
luminosity to reproduce this value. 

Abundances. The ratio of He to H was constrained from the 
relative strength of the optical He n lines (especially those at 
4686 and 5412 A) to the Balmer lines (mainly Up and U s ). 
The carbon abundance was obtained from C iv /L15801, 5812 
and, to a lesser extent, from C iv /Ml 548, 155 1 . The nitrogen 
content was derived from the same lines used to constrain 
the temperature. Note that the latter depends on the relative 
strength of those lines, while the N abundance affects the 
absolute strength of the lines. In addition, Niv /U486 and 
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Table 1. Photometry of the target stars from Massey & Duffy 
d2001l) together with the derived extinction. 



Star 


ST 


U 


B 


V 


E(B-V) 












Gal.+ SMC 


SMC-WR1 


WN3h 


14.22 


15.13 


15.17 


0.24 


SMC-WR2 


WN5h 


13.26 


14.08 


14.26 


0.03+0.06 


SMC-WR4 


WN6h 


12.48 


13.25 


13.37 


0.04+0.10 



N iv A1120 were also use as a secondary indicator. O v A5591 
and O v /15604 are used to provide upper limits on the oxygen 
content. 

• Mass loss rate and terminal velocity. M was derived from the 
numerous emission and P-Cygni lines present throughout the 
(F)UV and optical spectrum. We used the blue-ward velocity 
edge of the absorption part of P-Cygni profiles to derive Voo. 

• Velocity field. The slope of the velocity field (described by 
the so-called f3 parameter) was derived from the shape of the 
He n 4686 emission profile. This profile reacts to a combina- 
tion of change of density due to a different velocity field and 
to the subsequent change of formation depth of the line im- 
plied by the modified atmospheric density structure. In prac- 
tice, a narrower, more peaked profile results from a larger /3. 
This is somewhat different from the usual behavior encoun- 
tered among O stars. The difference is due to the higher wind 
density affecting the line formation depth in WR stars, while 
for O stars it is not significantly changed. 

• Clumping. The clumping factor /co was constrained from the 
strength o f the el e ctron scattering wing of Hen /14686. As 
shown by iHillierl ([1991), this part of the spectrum is espe- 
cially sensitive to the degree of homogeneity in the star's 
atmosphere: more structured winds lead to weake r electron 
scatter ing wings. This method was already used bv lCrowfherl 
d2000t) for the study of the star SMC-WR4. 

Finally, in the FUSE range we have added HI and H2 inter- 
stellar absorptions on top of our synthetic spectrum. We have 
not attempted a detailed fit of these interstellar lines. Instead, we 
have chosen column densities reproducing the main absorption 
lines. 

4. Results 

The results of our analysis in terms of derived parameters are 
summarized in Table|2] In the following, we briefly comment on 
the analysis of each individual star. 

4.1. SMC-WR4 

The fit of the SED is shown in Fig. Q] The error bars on the opti- 
cal data reflect the dispersion of measured V magnitudes found 
in the literature. A luminosity of 10 5 90 L Q is derived. The best fit 
of the normalized spectra of SMC-WR4 is shown in Fig. IA.1I 
The only major difficulty concerns the emission line red-ward 
of N v /14604, which is not reproduced. But overall, the fit is of 
excellent quality. An effective temperature of 42000 K gives the 
best fit of the optical N lines. Slightly lower values (~ 41000 
K) are preferred for the near-IR He lines (especially Hei A 1.08 
fim), still consistent with the N lines results within the uncer- 
tainties. This difference can be explained if the star's spectrum 
shows variability. Indeed, the near-IR and optical data have been 
obtained at different epochs. Since the star shows line profile 
variability (see paper I), simultaneous optical/infrared observa- 
tions are needed to investigate this discrepancy. The simplified 



approach to treat clumping might also affect the detailed line 
profiles. For the FUV interstellar absorption, we found that HI 
(resp. H2) column densities of 10 22 cm~ 2 (resp. 10 19 7 ) gave rea- 
sonable results. Individual abundances for He, C, N and O (up- 
per limit) were derived from our analysis. Their determination 
is illustrated in Figs. [2] A medium degree of He enrichment is 
detected in SMC-WR4, together with a rather strong N over- 
abundance and a carbon depletion. The typical uncertainties on 
these abundance determinations are of the order 30 to 50%, as 
seen from Fig. [20. Finally the mass loss rate determination gives 
a value of 8 x 10~ 6 M . Strong optical emission lines would re- 
quire slightly larger values of the order 10 M , while near-IR 
H and He lines tend to indicate slightly lower values ( ~ 7 x 10~ 6 

m ). a 

The left panel of Fig[3] shows various models with different 
clumping parameters (/«, ranging from 0.01 to 0.5). For each 
model, the mass-loss rate is scaled to properly account for the 
emission of He 11 /14686. One clearly sees that the strength of the 
red electron scattering wing of the line depends on the amount 
of structure in the wind. From this feature, we see that a model 
with fco ~ 0.1 best represents the observed profile. For a more 
quantitative statement, we have estimated the goodness of the fit 
of this red wing by means of a;^ 2 analysis of the region between 
4702 and 4730 A. The results are shown in Fig. [3] right panel. A 
value of 0.15+0.05 is preferred. 

SMC-WR4 was studied by ICrowtheJ d2000h in the optical 
and near infrared, using CMFGEN. The parameters we obtain 
are very close to his. The only notable difference is the slightly 
larger luminosity we derive, although within the error bars the 
agreement is acceptable. The difference comes from the differ- 
ent di stances to the SM C that were adopted: we use 63.2 kpc 
while Crowther (2000) assumes a distance of 60.2 kpc. Using 
his distance, we would derive a luminosity of 10 5 8 L Q , in good 
agreement with Crowther's result. We also note that Crowther's 
fit also showed the problem reported for the mass loss rate de- 
termination: near-IR lines tend to indicate lower M than optical 
lines. 

4.2. SMC-WR2 

The fit of the SED of SMC-WR2 is shown in Fig. Q] The er- 
ror bars on the optical data reflect the dispersion of measured V 
magnitudes. We derive a luminosity of 1O 5 5O L . Fig. lA.2l shows 
our best fit of the normalized spectra. The overall quality of the 
fit is as good as for SMC-WR4. HI and H2 interstellar absorp- 
tion lines have been added in the FUSE range with the follow- 
ing column densities: 10 2L7 cw -2 for HI and 10 193 cm~ 2 for H2. 
The optical N lines indicate an effective temperature of ~ 42500 
K. We note that with this T e fr the Hei /1 1.08 fim line is under- 
predicted. Obtaining a stronger emission requires a reduction 
of T e jf by ~ 5000 K. But in that case the optical spectrum is 
not reproduced anymore: the Nv lines vanish completely, Niv 
/14058 A is strongly reduced and the Nm lines become much 
too strong. The same explanation as for SMC-WR4 could be 
invoked. Clearly, coordinated multiwavelength observations (as 
soon possible with instruments such as VLT/X-Shooter) are nec- 
essary to solve this issue. One can also note that the S/N ratio of 



1 These errors do not include any systematics due to uncertainties in 
the atomic data 

2 Note that our best fit model in Fig. I A. ll has a slightly smaller mass 
loss rate (0. 1 dex) than the best fit model of Fig.[2]in order to best match 
both optical/UV and near-infrared lines. In Fig.(2] M is scaled to fit He 11 
,14686 
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2000 3000 4000 5000 6000 2000 3000 4000 5000 



A[A] A[A] 

Fig. 1. SED fitting for SMC-WR2 (left panel) and SMC-WR4 (right panel). The IUE spectra are reported as well as the fluxes in 
the UBV bands. E(B-V)=0.03+ 0.08 (Galactic + SM C) and 0.04+0.10 are required for SMC-WR2 and SMC-WR4 respectively. A 
distance of 63.2kpc is assumed ddi Benedettoll 1997b . This SED fitting yields the luminosities: 1O 5 - 5O L for SMC-WR2, 1O 5 9O L for 
SMC-WR4. 




Fig. 2. He (left) and CN (right) abundance determination for SMC-WR4. The black solid line is the observed spectrum, the inter- 
rupted colored lines are models with different He/H (C/H, N/H) ratios. 



the near-IR data is lower than for the optical data, increasing the 
uncertainty on the temperature determination from NIR lines. 
Abundance determinations were performed as for SMC-WR4. 
Only a small He enrichment is required to fit the He and H lines 
of the spectrum. However, evidence for CNO processing comes 
from the existence of a C depletion and a N overabundance. We 
note that both optical and near-IR emission lines are well repro- 
duced by the a single value of the mass loss rate. 

Unlike in SMC-WR4, He n /14686 is a rather weak emission 
line. Consequently, its red wing is more sensitive to details of the 
data reduction and normalization. We have attempted to derive 
the clumping factor /„, as for SMC-WR4, but it turned out that 
no conclusions could be drawn due to these observational un- 
certainties. In particular, the exact position of the continuum is 
uncertain due to residuals of the flat field correction and to im- 
perfect order merging. Other classical clumping indicators are 



found in the FUV-UV range. Ov /1 1371 and Niv ^1720 have 
been shown to dep end on the degree of inhomogeneities (e.g. 
iBouret et al.1 120051) . However, the IUE spectra of SMC-WR2 
have a too low spectral resolution and these lines cannot be used. 
In the FUSE range, P v AA1 1 18, 1 128 are kno wn to depend on 
the amount of clumping. Full erton et alj d2006l) showed that very 
low filling factors we re necessary to accoun t for the strength 
of this doublet, while ICrowther et alJ (l2002al) a nd iHillier et al. 
(2003) obtained satisfactory fits with /«, ~ 0.1. Fullerton et al. 
d2006 ) highlighted the high sensitivity of this doublet to the 
ionization structure, which can be influenced by shocks and 
X-rays. This effect is probably reduced in dense WR star 
winds compared to O stars, but might still be present. Besides, 
Oski nova et alJ (120071) revealed that depending on which type 
of clumping was present (micro vs m acro-clumping) the line 
profiles had different shapes. Finally, ICrowther et alJ ((2002a) 
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Fig. 3. Determination of clumping parameter / M . Left: Effect of 1 
4686 for star SMC-WR4. Right: Quantification of the goodness 
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clumping parameter f ro on the electron scattering wing of Hell 
the fit of the red wing of Hell 4686. The solid line is a simple 
clumping factor / ro . 



showed that the phosphorus abundance strongly influenced the 
shape of this unsaturated line. Given these multiple dependen- 
cies, we refrained from using Pv /L11118, 1128 as a clumping 
diagnostic in the present study. We simply adopted /oo=0.1 for 
SMC-WR2. 



4.3. SMC-WR1 

Unfortunately, only optical and near-IR spectra are available for 
SMC-WR1. The analysis is thus more limited than for the two 
other stars. In particular, in the absence of UV P-Cygni lines, we 
adopted a terminal velocity of 1800 km s _1 which turned out to 
give an excellent fit of the optical emission lines. The uncertainty 
on T e ff (and consequently on L) is larger for this star since only a 
limited number of lines from successive ionization states is avail- 
able. In particular, although several N v lines are used, the only 
N iv line (at 4058 A) is very weak. It can be used to place a lower 
limit on T e ff. The upper limit comes from the shape of the N v 
AA46Q4, 4620 lines which become narrower and weaker at too 
high effective temperature. In practice, the uncertainty on T e ft is 
+ 7000 K. This translates into an uncertainty on log -f- of 0.15 
dex. In spite of these increased uncertainties, we were able to 
derive a number of stellar and wind parameters. Fig. IA. 31 shows 
our best fit. Overall, the quality is good: optical and near-IR lines 
can be reproduced with a single set of stellar and wind param- 
eters. The He n lines of the Brackett series are slightly too deep 
and narrow. A change of f3 does not help since it degrades the 
fit of Hen A4686. Similarly, increasing v M broadens too much 
the other emission lines. Hen /15412 is slightly underestimated. 
Improving its fit by increasing M, T e ff or He/H alters the fits of 
all other lines. Presumably, remaining uncertainties in the den- 
sity structure are responsible for this discrepancy. 

Like SMC-WR2 and SMC-WR4, SMC-WR1 is slightly He 
enriched and strongly N rich. Fig.|4]illustrates our He/H determi- 
nation. As we will see later, this quantity is especially important 
to understand the evolution of SMC-WR1 . From the weakness of 
C iv and O v lines in the red optical spectrum, one can put upper 
limits on the C and O content (see Table |2), indicating deple- 
tion of these two elements. As for SMC-WR2, He n A4686 is not 
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Fig. 4. Determination of the He to H ratio for star SMC-WR1. 
The observed spectrum is the black solid line. The mass loss 
rates of the models with different He/H ratios were adjusted to 
give a correct fit of the He n /14686 A emission line. A value of 
He/H=0.25 is preferred. 



strong enough to allow a correct clumping factor determination. 
We adopted =0.1. 



5. Discussion 

In the next sections, we first comment on our results on the 
clumping factor of SMC-WR4 in the context of a metallicity de- 
pendence of clumping properties. We then move to the interpre- 
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Table 2. Parameters of our best fit models. Uncertainties are: 3kK (7kK for SMC-WR1) on temperatures, 0.1 dex on luminosities 
(0.15 for SMC-WR1), 0.2 dex on M, 100 km s 1 on v M , 0.05 on He/H, 30 to 50% on abundances. 



Star 


ST 


T, 


T c ff 




R, 


logM 


f« 


P 


v„ 


He/H 


X 


X(C) 


X(N) 


X(O) 






[kK] 


[kK] 


Ro 


[Mo yr" 1 ] 






[km s" 1 ] 


# 


(%) 


(%) 


(%) 


(%) 


SMC-WR1 


WN3h 


65.4 


60.5 


5.75 


5.8 


-5.55 


0.1(£ 


1.0 


1800 


0.25 


0.50 


<0.004 


0.12 


<0.04 


SMC-WR2 


WN5h 


44.8 


42.5 


5.50 


9.3 


-5.85 


0.10~ 


1.0 


1000 


0.15 


0.62 


0.005 


0.35 




SMC-WR4 


WN6h 


43.1 


42.0 


5.90 


16.0 


-5.00 


0.15 6 


1.0 


1000 


0.45 


0.35 


0.003 


0.40 


<0.15 



" Adopted 

* Derived from^ 2 fit, mass loss rate scaled from model closest to minimum (^=0.1). 



tation of the evolutionary status of the three WNh stars analyzed 
in the present study. 

5.1. Winds and clumping 

In paper I, a thorough investigation of moving features in strong 
emission lines revealed that the dynamical properties of clumps 
in the SMC were very similar to those of Galactic WR stars. 
Here, we have attempted to constrain the clumping factor at low 
Z. For star SMC-WR4, we have carried out a detailed study of 
the clumping factor. We have quantitatively shown that the pre- 
ferred value for f x (see Eq. [TJ was 0.15+0.05. The^ 2 curve of 
Fig.[3]unambiguously confirms that lower and larger values are 
excluded. 

Studies of Galactic WR stars have been the first to re- 
veal that clumping was present in m assive-star winds (e.g. 
iMoffat et al.lll988t Ist-Louis et al.lll989b . Quantitative determi- 
nations of clumping factors were carried out by various groups 
using different atmosphere codes. lHamann & Koesterke ( 1998) 
analyzed three WN (one WN5, one WN6 and one WN8) stars 
and one WC star and found that a clumping filling-factor of 0.25 
was best suited for WN stars, although 0.06 was still compat- 
ible within the unc ertainties. The latter value was best for the 
WC star. Similarly, ISchmutzl d 19971) found that /«, = 0.25 led 
to a good fit of the em ission lines of the WN4 star HD50896. 
iHillier & Miliel d 19991) found that / M = 0. 1 was bes t suited to fit 
the UV -optical-IR spectrum of a Galactic WC5 star. lHillier et alj 
(2001) carried out an extensive study of the spectra of 77 Car 
and found that a clumping factor of ~0. 1 better fitted the elec- 
tron sratterin^wings_of both H and Hel lines compared to ~0.2- 
0.3. lHerald et all d2001l) analyzed two Galactic WN8 stars and 
concluded that 0.05 < / < 0.2 from the same diagnostics. 
ICrowther etafl d2002al 120061) analyzed WC stars in the LMC 
and the Galaxy and found that f ra =0. 1 was preferred compared to 
1 . and . 1 . Howev er, no other intermediate values were tested. 
Naia rro et"aT] d2004l) found values between 0.08 and 0. 15 for 
five WNLh stars in the Arches cluster. Hence, a value of 0. 1-0.2 
seems to be common among Galactic WR stars and related, post 
main sequence objects. 

For SMC-WR4, we have shown that f x was definitely be- 
tween 0.10 and 0.20 (see FigOJl. Keeping in mind that we are 
dealing with a single measurement for the SMC, it looks as if 
there is no difference between the low metallicity environment 
of the SMC and the Galaxy in terms of clumping factor. Or more 
accurately, if there is any variation, it is not larger than ~ 0.1, 
which is the uncertainty of current determinations. This is con- 
sistent with the results of paper I. 

From a theoretical point o f view, the hydrodynamical s imu- 
lations of[O wocki et al. ( 1988), R unacres & Owockl d2002l) and 
Runacres & Owockil d2005h show that instabilities grow at a rate 
proportional to g m d, the line radiative acceleration. The radia- 



tive acceleration g mc i scales with Z 1 ~" (see Eq. 12 and 79 of 
iPuls et al]|2000l) . a being the classical CAK parameter. But g ra d 
also depends on p ~ a , and since p (the density) is directly pro- 
porti onal to M, wh ich in turn scales roughly as (Eq. 87 
of IPuls et al.l l2000), all Z effects cancel out. In this context, no 
metallicity effects are expected if clumping is driven by rad iative 
linear instabilities. Alternatively, Gay l ey & Owockil dl995) have 
shown that in the case of multiple scattering, line overlap leads 
to a reduction of the linear instability growth rate. Presumably, 
lowering the metallicity would reduce this overlap and produce 
an enhanced growth rate (S. Owocki, priv. comm.). In that case, 
clumping might be stronger at low Z. This is however still widely 
speculative, and future hydrodynamical simulations, including 
non linear effects, will certainly shed more light on this issue. 

5.2. Evolutionary status 

We have shown that the three WNh stars analyzed are slightly He 
rich, strongly N rich, and C (and O for SMC-WR1 and SMC- 
WR4) poor. This qualitatively indicates that they are massive 
stars in an early phase of their post-MS evolution. The large 
quantity of H still present in their atmosphere shows that they 
are still H-burning objects, or that they have entered the core 
He-burning phase very recently. This is confirmed by the posi- 
tion of the stars in a diagram showing the H mass fraction as 
a function of luminosity (Fig. [5] left): one sees that the three 
stars are located on the early parts of evolutionary tracks, still 
far away from classical H-free WR stars. In this diagram, all 
stars can in principle be accounted for by evolutionary tracks. 
This was not th e case in the earlier study of SMC-WR4 by 
ICrowfherj (;2000): he found that SMC-WR4 was underluminous 
by rou ghly 50% compared to the closest track of M evnet et alj 
dl994l) . Crowther concluded that stellar evolutionary tracks were 
probably missing an important ingredient. This ingredient was 
stellar rotation, wh i ch is included in the evolutionary tracks of 
Mev net & Mae der (2005) used in Fig. (the initial rotational 
velocity in these tracks is 300 km s ). From this diagram, one 
can estimate the initial masses of our program stars by interpo- 
lating between existing tracks (at Z = 0.2Z o ). The results are 
given in the first column of Table[3] masses range from about 25 
M Q (SMC-WR2) to about 40 M (SMC-WR4, SMC-WR1). The 
errors take into account the uncertainties on T e ff and log jj-. 

The middle panel of Fig. [5] shows the C versus N ma s s frac- 
tion together with the predictions of Meynet & Maeder (2005) 
for various initial masses. The errors on the carbon content are 
smaller than the symbols. SMC-WR4 and SMC-WR2 are in (or 
very close to) a phase where X(C) is minimum and X(N) is max- 
imum. SMC-WR1 appears to have a carbon mass fraction con- 
sistent with CNO equilibrium, but its N content is a little too 
low. A similar conclusion is reached if one looks at the right 
panel of Fig. |5j it shows the O versus N mass fraction for stars 
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Table 3. Initial mass estimate (in solar units) from evolutionary 
diagrams. The second and third column give estimates based re- 
spectively on the X(H)-L and HR diagram. The errors take into 
account uncertainties on T e ff and log j-. Tracks with Z = 0.2Z o 
from lMevnet & Maederl d2005l) are used. 



Star X(H)-L T cff -L 



1 


34.4^' 


>40 


2 


22.8+^ 


47 3 +7 -° 


4 


4.7 7+6.9 
^"-7J 


64 7 +16 ° 



SMC-WR1 and SMC-WR4. As for carbon, the oxygen content 
is reduced during the CNO cycle, so that an O underabundance 
can be expected from evolutionary tracks. The O mass fraction 
of SMC-WR1 is consistent with the expected minimum, as is the 
case for carbon. In this figure too the N mass fraction is not as 
large as expected for the O content. We will see later that the 
current tracks might not be appropriate for this star. The upper 
limit on the O content of SMC-WR4 is too high to draw any 
meaningful conclusion. This limit (not shown) is even higher for 
star SMC-WR2. 

Let us now turn to the HR diagrams shown in Fig. [6] The 
position of the three ta rget stars is indicated on top of the 
iMevnet & Maeder (2005) evolutionary tracks at Z-0.2 Z (ap- 
propriate for the SMC, left panel) and Z=0.4 Z (right panel). 
We have highlighted in bold font the phase of the tracks in which 
the H mass fraction is larger than 20%. Since SMC-WR1, SMC- 
WR2 and SMC-WR4 all have X(H) > 0.3, any comparison be- 
tween the position of the stars and evolutionary models should 
rely on the bold lines in either panel. Let us consider each star 
one by one: 

- SMC-WR4: for this star, we find X(H)=0.35±0.1. From the 
left panel of Fig. [6] (i.e. the one with evolutionary tracks at 
0.2 Z ), one might invoke two scenarios to explain the po- 
sition of SMC-WR4. In the first one SMC-WR4 is a rela- 
tively unevolved star, still close to the main sequence and 
with an initial mass slightly above 60 M . In that case the 
initial mass can be estimated from the uncertainties on T e ff 
and log y-, as given in the second column of Table [3] The 
result is somewhat surprising, since it is much larger than 
the estimate based on the X(H)-L diagram, with no overlap 
within the errors (42.2+^ vs 64.7^°). Besides, if SMC- 
WR4 were indeed at an early post-main sequence phase of 
its evolution, its mass loss rate wo uld be around 10~ 5 7 M 
yr~' according to lVink et al.l (12001). Instead, we find a value 
one order of magnitude larger (once corrected for clumping). 
Thus it appears unlikely that SMC-WR4 is just evolving off 
the main sequence. 

Another possibility is that SMC-WR4 is a star with an initial 
mass of 40-50 M and evolving blue-ward from a supergiant 
phase. Indeed, the Z = O.2Z 40 M track ends a little before 
reaching the position of SMC-WR4. However, the same 
track at Z=0.4 Z passes through the position of SMC-WR4 
while still having X(H)>0.2 (see right panel of Fig. [SJ. 
Hence assuming a metal content slightly larger than 0.2 Z 
can explain the characteristics of SMC-WR4. In that case, 
the 40 M track would explain the position of SMC-WR4 
in both the HR and the X(H)-L diagram (the initial mass 



derived from the latter diagram is very similar at 0.2 and 
0.4 Z for SMC-WR4). Note that other effects such as a 
slightly different initial rotation rate com pared to the values 
used i n the evolutionary calculations of M evnet & Maederl 
(120031) . or different mass loss rates (uncertainties being large 
due to clumping) could also modify the tracks at Z=0.2 
Z and explain the position of SMC-WR4. Given this, we 
tentatively conclude that SMC-WR4 results from rather 
normal evolution. 

- SMC-WR2: the hydrogen mass fraction is the largest of the 
three stars, reaching 0.62+ 0.1. The situation is different 
from SMC-WR4. There is only one possibility to explain 
both the effective temperature/luminosity and H content 
of SMC-WR2: it must be a star just evolving off the main 
sequence. Indeed, none of the evolutionary tracks coming 
back from the cool part of the HR diagram and with low 
enough luminosities can account for the large H mass 
fraction of SMC-WR2. The 3OM track at Z=O.4Z has a 
too low H content when reaching the position of SMC-WR2 
(see the tick on the track at X(H)=0.3). The star has to be 
in its early evolution. If we accept this, one can estimate the 
initial mass of the star by interpolation between the red-ward 
evolving tracks. The results (given in the third column of 
Fig. [3]) indicate M=47.3+j ° 6 M . This is significantly larger 
than the estimate from the X(H)-L diagram (22.8^ M ). 
Such a difference is surprising, since it is based on the 
same set of tracks. We will discuss possible solutions to this 
puzzle in Sect. 15.31 

- SMC-WR1: The situation is even more intriguing for star 
SMC-WR1. We find convincing evidence that SMC-WR1 is 
still H-rich, with X(H)=0.50+0.1. But a look at Fig.|6]shows 
that SMC-WR1 is definitely located to the left of the H- 
burning main sequence, at a temperature of ~65000 K and 
log y- of 5.75. The only evolutionary tracks reaching this 
part of the diagram have initial masses larger than ~ 50 M . 
But this is only at the very end of the star's evolution, when 
it is H and N free. The properties of SMC-WR1 are thus ex- 
tremely disturbing, classical stellar evolution being unable to 
explain them. 



From this analysis, one can conclude that one star (SMC- 
WR4) seems to have followed a normal evolution. On the con- 
trary, the two remaining stars (SMC-WR1 and SMC-WR2) have 
probably been affected by uncommon physical processes. We 
tentatively identify some of them in Sect. 15.31 

Before concluding this section, we briefly note that mass es- 
timates are important clues to the nature of the sample stars. 
One can add another such estimate assuming that the stars have 
Eddington luminosities. In that case one obtains lower limits to 
the current masses of the three WNh stars. These values are re- 
spectively 9.8, 7.8 and 19.6 M for SMC-WR1, SMC-WR2 and 
SMC-WR4 (assuming the only source of opacities is electron 
scattering). In practice, our modelling indicate that the maxi- 
mum Eddington factors in their atmosphere are of the order 0.7 
and that the total radiative acceleration is 1.2-1.5 times the ac- 
celeration due to electron scattering, implying current masses of 
the order 10-35 M . Of course, these values have to be taken as 
rough estimates since they depend on the exact hydrodynamical 
structure which still suffers from uncertainties (see Sect. [3). 
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Fig. 5. Evolutionary status from chemical composition. Left: Hydrogen mass fraction as a function of Luminosity. Middle: Carbon 
mass fraction as a function of Nitrogen mass fraction for stars SMC-WR2 and SMC-WR4. Right: O mass fraction as a function of N 
mass fraction from evolutionary tracks (lines) and as derive d for star SMC-WR1. The three stars are still rather H rich, but already 
show signs of CNO processing. Evolutionary tracks are from lMevnet & Maederl (120051) and have an initial rotational velocity of 300 
km s . 





Fig. 6. HR diagram with the three sample stars. Evolutionary tracks are from iMevne t & Maeder (2005), the left panel being for 
Z=O.2Z , the right one for Z=O.4Z . We have highlighted the parts of the tracks corresponding to the H-rich phase, i.e. X(H)>20%, 
using bold lines. On the Z=O.4Z 30 M Q track, we have also noted the position where X(H)=0.3. 



5.3. Evidence for evolution with fast rotation. 

The previous section has shown that classical stellar evolution 
was able to reproduce the observed properties of SMC-WR4 if 
its metallicity is a little larger than 0.2Z o . However, it was shown 
that the current set of evolutionary tracks failed to explain SMC- 
WR1 and SMC-WR2. 

Before considering SMC-WR1 and SMC-WR2, let us first 
reiterate our conclusions regarding SMC-WR4. We showed that 
both its position in the HR diagram and its surface abundances 
were compatible with a star of initial mass 40-50 M Q in or about 
to enter the core He burning phase. This is rather similar to 
the conclusion of iFoellmi et al.1 J2003) who argued that SMC- 
WR4 most likely resulted from standard evolution of a 70-85 
M s star. We stress once more that according to this scenario, 
SMC-WR4 is coming back from the red part of the HR diagram. 
Consequently, it is not a core H burning object. This is important 



in the context of the understanding of the properties of H rich 
WN stars. Recently, iMartins et aTJ(l2008l) analyzed 16 WN7- 
9h stars in the Arches cluster and showed that they were rather 
young and very massive stars still in the core H burning phase. 
The WN 6h stars in the core of NG C3603 are probably similar 
objects dCrowther & Dessartlll998l) . The initial mass estimates 
for these objects in the A rches and NGC3603 cluster range be- 
tween 60 and 130 M Q ("see lSchnurr et al.l(l2008l) for a recent kep- 
lerian mass estimate of NGC3603/A1). Hence, it seems that H- 
rich WN stars come in different flavors. They can be very mas- 
sive objects in an early phase of evolution, with ages not larger 
than 2-3 Myr; but they can also result from less massive objects 
and be more evolved. In the latter category, one should men- 
tion st ar WR3 in the outer part of the Galaxy. Marchenko et al.l 
(2004) showed that it was a H-rich object resulting from the evo- 
lution of an initially 40-50 M object. Since lower metallicity is 
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encountered in the outer regions of the Galaxy, star WR3 might 
be more similar to the SMC stars analyzed in the present study. 

Coming back to SMC-WR1 and SMC-WR2, there are a few 
possibilities to e xplain the puzzling r e sults o f Sect. 15.21 The first 
one is binarity. IWellstein & Langerl ([T999) showed that mass 
transfer during binary evolution could dramatically change the 
evolutionary sequence of both components of the system. In par- 
ticular, their Fig. 1 shows that the primary might evolve almost 
vertically just after the mass transfer. Consequently, if such a 
star was analyzed by means of classical single star tracks, its 
mass would be overestimated. This could explain the inconsis- 
tencies encountered for SMC-WR2: its mass derived from the 
HR diagram is larger than that derived from the X(H)-L di- 
agram. Besides, mass transfer would affect the surface abun- 
dances, so that the derived H mass fraction might be different 
from that of single evolutionary sequences. Hence, mass esti- 
mates from the X(H)-L diagram would be affected too. Binarity 
might thus be a n alternative to explain the properties of SMC- 
WR2. However, Foell miet al.l d2003l) thoroughly studied the ra- 
dial velocities and photometry of all the SMC stars and showed 
that SMC-WR2 was not varying. They concluded that it was 
most likely a single star, or maybe a long period binary. But 
in any case it was not a short period binary prone to experi- 
ence the kind of evolution described above. The same conclu- 
sion was reached for SMC-WR1. Periodic changes of the shape 
of emission profiles may also reveal the close-binary nature of 
a massive star. The line-profile variability of SMC-WR1 and 
SMC-WR2 closely resembles the behavior of emission lines in 
the presumably single Galactic W-R stars (Paper I), where the 
stochastic, seemingly random variability is induced by numerous 
wind-embedded clumps. Hence, we conclude that practically all 
the variability patterns seen in SMC-WR1 and SMC-WR2 are 
caused by stochastic appearance of small-scale overdensities in 
the winds of the SMC W-R stars. 

Binarity thus appears unlikely to be the explanation for the 
properties of SMC-WR1 and SMC-WR2. Alternatively, there 
is growing evidence that homogeneous evolution might be at 
work in at least a fraction of massive stars in the Magellanic 
Clouds. IBouret et al.l d2003l) argued that star MPG 355 in NGC 
346 could have the same age a s the other clu s ter me mbers if it 
was evolving homogeneously. Walbo nTet al.l (|2004) explained 
the position of highly N enriched Magellanic Cloud 02 stars in 
the HR diagr am by large mixing le ading to homogeneous evolu- 
tion. Finally. iMokiem et al.l d2007al) explained the correlation of 
large mass discrepancies Q with He enrichment by the fact that 
He-rich stars might evolve homogeneously, so that using normal 
tracks would overestimate their evolutionary mass. 

Homogeneous evolution w as shown to occur in ra pidly rotat- 
ing sta rs half a centu ry ago by ISchwarzschildl £l958). More re- 
cently, [MaedeB (EUlb and lLangeridl992l) ran evolutionary mod- 
els confirming that fully mixed stars were evolving blue-ward 
in the HR diagram. The conditions under which such homoge- 
neous chemical evolution occurs depend on several parameters. 
The basic constraint is that the timescale for rotationally induced 
mixing should be shorter than the nuclear timescale. In practice, 
this means that the stars should keep a high rotational rate dur- 
ing their evolution. A key ingredient to keep the high rotational 
rate is the strength of the stellar wind. Larger mass loss rates 
will generally lead to larger angular momentum removal, which 



3 The "mass discrepancy" problem refers to the larger evolutionary 
masses derived from HR diagrams compared to spectro scopic masses 
derive d from gravities through atmosphere modelling (see Herre ro et alj 
1992). 



in turn will slow down the rotation. As radiatively dri ven winds 
are weaker at lower metallicity (e.g. iPuls et af]|2000l) . changes 
in the surface abundances are expected to be more common 
at (very) low metallicity. But the downward re vision of mass 
loss ra t es of massive stars due to clumping (e.g. ICrowther et all 
2002ai IBouret et al.ll2005l) could make homogeneous evolution 
more common at higher metallicity than previously thought. 
Besides, mass loss anisotropy implied by fas t rotation might af- 
fect t he details of the subsequent evolution (Meynet & Maeder 
120071) . In addition, magnetic fields can modify the way rotation 
chang es during the lifetime of t he star, leading to solid body ro- 
tation dMevnet & Maederll2003l) . From this, one sees that homo- 
geneous chemical evolution will occur under a combination of 
effects. A systematic study of these conditions is still lacking, 
partly because of the uncertainties in the mass loss rates of mas- 
sive stars, especially Wolf-Rayet stars. 

Nonetheless, various works have mentioned the existence 
of homogeneous evolution. Meynet & Maeder ( 2000) argue that 
stars initially more massive than 40 M Q and with > 0.5 - 
Q (Qen,) being the (critical) rotation rate - will evolve homo- 
geneo usly. The 60 M evolutionary track of lMeyneT & Maeder 
(120051) with an initi al V sin; of 500 km s" 1 , as well as the fast 
rotating models of Meynet & Maeder (2007) (w hich include the 
effects of a magnetic field) confirm this trend. lYoon & L anger 
(2005) also produce blue-ward evolving stars in their fast rotat- 
ing models at low metallicity. As mentioned above, a complete 
set of evolutionary tracks for various metallicities and rotational 
rates does not exist at present. Only a few tracks are available. 
In Fig. [7J we plot t he ev olutionary tracks of a 60 M Q star from 
lMevnet & Maeder! d2005l) for both V sin/ = 300 km s" 1 (solid 
line) and 500 km s (dashed line). The latter corresponds to a 
borderline case where the evolution is almost homogeneous, the 
star evolving at almost constant T e ff in the first phases. In the 
X(H)-L diagram, one sees that for fast rotating stars, the track 
is shifted to lower luminosities. In practice, this means that the 
mass one would derived from a set of tracks with V sin/ = 500 
km s would be larger than for the case where V sin/ = 300 km 
s _1 . On the contrary, fast rotating stars are on average more lumi- 
nous in the HR diagram (due to the effects summarized above). 
Hence, lower masses are derived for larger V sin/ tracks. This 
is interesting in the context of the mass problem reported for 
star SMC-WR2: the large difference between the masses derived 
from the X(H)-L and HR diagram can, at least qualitatively, be 
reduced if one uses faster rotating tracks. This strengthens the 
idea that SMC-WR2 might be a rapidly rotating star, close to 
homogeneous evolution. It is worth mentioning that using the 
two methods to estimate masses we have applied in this paper, 
one might have an indirect way of determining the initial rota- 
tional velocity of Wolf-Rayet stars: a set of rotating tracks should 
be able to give consistent mass estimates in the X(H)-L and HR 
diagram. If not, the initial rotational velocity is probably not cor- 
rect. 

The case of SMC-WR1 is more extreme: it is located on the 
left side of the main sequence, but still co ntains a significant 
fractio n of hydrogen. Inspection of Fig. 1 of iMevnet & Maeden 
(2007) reveals that homogeneous evolution can reproduce these 
features. In this figure, SMC-WR1 would lie close to to the 
anisotropic wind track (log T e fi = 4.8, log ^- = 5.75) - although 
slightly below due to the lower mass of the star compared to 
the model -, where the core (and consequently the surface) still 
contains hydrogen. Qualitatively, SMC-WR1 is thus likely an 
example of a WR star proceeding from homogeneous evolution 
of a massive star. Star WR3 in the Galaxy is probably a twin 
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Fig. 7. Illustration of the effect of rotation on the X(H)-L diagram (left panel) and on the HR diagra m (right panel). The masse s 
estimated using different rotational velocities will be different in both diagrams. Evolutionary tracks from Mevnet & Maeder (2005). 



of SMC-WR1 . iMarchenko et all (12004 showed that it was H- 
rich and very similar spectroscopically to early WN stars in the 
SMC. Given our results, it is likely that WR3 also results from 
evolution with high rotation. 

As stated above, no systematic study of the occurrence of ho- 
mogeneous chemical evolution with metallicity, magnetic field 
and rotation velocity exists. However, a key ingredient is the 
high initial rotation rate. One might wonder if initial rotational 
velocities as large as 500 km s _1 can be reached in the SMC. 
Recent studies of the projected rot ational velocitie s of B stars 
in the SMC were carried out by Marta van et al.l d2007l) and 
iHunter et al.l d2008l) . Both works show that the V sin; distribu- 
tion is shifted towards larger values in the SMC compared to 
th e Galaxy, with ave rages around 155-175 km s _1 . According 
to lHunter etail (120081) . about 10% of the SMC stars have V sin; 
larger than 300 km s _1 . Considering that these values are lower 
limits on the true surface velocities (the inclinations being un- 
known), stars rotating at 500 km s" 1 might exist in the SMC. We 
thus argue that SMC-WR1 is most likely the result of homoge- 
neous evolution of a massive star. 

We have seen that one could get insight into the initial ro- 
tational velocity of our target stars by using both the X(H)- 
L and HR diagrams. In Fig. [8] we constrain the current pro- 
jected rotational velocities of SMC-WR1 and SMC-WR2 from 
the shape of narrow absorption/emission lines. One clearly sees 
that values much larger than 50 km s _1 are excluded, the line 
profiles becoming much too broad. These values are lower lim- 
its on the true rotational velocities. However, it is unlikely 
that we see both stars pole-on, so that SMC-WR1 and SMC- 
WR2 probably have equatorial velocities < 100 km s -1 . Normal 
evolutionary tracks with rotation predict equatorial velocities 
betwe en and 200 kin s" 1 f or 40-60 M stars in the WR 
phase (iMevnet & Maeder 2005). On the ot her hand, the homoge- 
neous tracks of Mevnet & Maeder) (120071) predict velocities be- 
tween 100 and 400 km s in the case of isotropic winds, and 
around 600 km s _1 for anisotropic winds. The latter predictions 
(anisotropic winds) are supposed to be more realistic. Since we 
observe rather modest rotational velocities, and given the strong 
evidence for homogeneous evolution - and thus fast initial rota- 



tion - at least for star SMC-WR1, one might argue that braking 
was more efficient than predicted for these stars. 

In the con text of GRB formation, f ast rotation is the key in- 
gredient (see lWooslev & Bloom! (l2006h for a review). The other 
important condition is that the star explodes while it has lost 
all its hydrogen envelope. Here, we have shown that SMC-WR1 
has most likely evolved with a high rotation rate until recently. 
The fact that it still contains a large H mass fraction might ar- 
gue against its potential explos ion as a GRB. However, a ccord - 
ing to the model predictions of M evnet & Maedej d2007l) . a star 
following homogeneous evolution will end its lif e as a H-free 
WR star (of the class WC or WO). iHirschi et all (120051) argue 
that GRBs can exist at SMC metalliciti es and that t hey sh ould 
have WO stars as progenitors. Similarly. I Yoon et all J2006) pre- 
dict a metallicity threshold of about Z=0.004 for the formation 
of long-soft GRBs (the exact value depending mainly on the 
adopted mass loss rates). In their scenario, stars ending then- 
lives as GRBs follow quasi chemically homogeneous evolution. 
The properties of SMC-WR1 are very reminiscent of this sce- 
nario. However, the questions raised in the the previous para- 
graph prevent us from concluding firmly that SMC-WR1 is on 
its way to become a GRB progenitor. 



6. Conclusion 

We have conducted a detailed spectroscopic analysis of three 
WNh stars in the SMC. Using atmosphere models computed 
with the code CMFGEN, we have determined the stellar and 
wind parameters. The key results can be summarized as follows: 

• the quantitative analysis confirms the spectral classification 
in the sense that all three stars still contain large amounts of 
hydrogen in their atmosphere. At the same time, they show 
clear signs of CNO processing, with C (and O for SMC- 
WR1) depletion, and N enrichment. SMC-WR1 is a hotter 
object than the SMC-WR2 and SMC-WR4, as suggested by 
its earlier spectral type. 

• SMC-WR4 can be explained by normal evolutionary tracks 
including rotation if its global metallicity is slightly above 
1/5''' solar. In that case, it is most likely a star with an initial 
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Fig. 8. Estimate of the current projected rotational velocity of 
stars SMC-WR1 and SMC-WR2. The black solid lines are the 
observed spectra of SMC-WR2 (top panels) and SMC-WR1 
(bottom panel). Colored lines are synthetic spectra convolved 
with rotational profiles with different V sin/. Values larger than 
50 km s are excluded. 

mass of 40 - 50 M in or very close to the core He burning 
phase. 

• SMC-WR2 and, more dramatically, SMC-WR1 cannot be 
accounted for by moderately rotating (i.e. 300 km s _1 ) evo- 
lutionary tracks. Their position in both the HR and X(H)-L 
diagrams can be reproduced only if the stars rotate initially 
very fast (> 500 km s _1 ). In the case of SMC-WR1, only ho- 
mogeneous evolution can account for the position of the star 
to the left of the main sequence and the high H mass fraction 
in its atmosphere. 

• the clumping factor /„, of star SMC-WR4 is quantitatively 
constrained to be 0.15±0.05. Within this uncertainty, no dif- 
ference is found with Galactic WR stars, confirming the dy- 
namical results of paper I. 

The present results are based on a limited number of stars. 
In a subsequent publication, we will present a study of the re- 
maining Wolf-Rayet population of the SMC. This will provide 
a wider view of the stellar properties of these objects. In partic- 
ular, we will be able to estimate the fraction of stars potentially 
following evolution governed by fast rotation. This will be use- 
ful in the context of Gamma-Ray Burst formation, thought to be 
due to very rapidly rotating WR stars. 

A more systematic study of the clumping properties of SMC 
WNh stars is also required to confirm that clumping does not 
seem to depend on metallicity. This is crucial to understand the 
physics of inhomogeneity formation in massive-stars winds and, 
more generally, the mechanisms of mass loss. Our future study 
will tackle these questions too. 

Acknowledgements. We thank the referee, Paul Crowther, for useful comments 
which helped to improve the spectroscopic analysis and for sharing near-infrared 
spectra. We thank Stan Owocki for illuminating discussions about radiative insta- 
bilities. FM thanks the CINES for allocation of computing time. AFJM thanks 



NSERC (Canada) and FQRNT (Quebec) for financial aid. JCB acknowledges 
support from the French National Research Agency (ANR) through program 
number ANR-06-BLAN-0105. 



References 

Blomme, R., Prinja, R. K., Runacres, M. C, & Colley, S. 2002, A&A, 382, 921 
Bouchet, P., Lequeux, J., Maurice, E., Prevot, L., & Prevot-Burnichon, M. L. 

1985, A&A, 149, 330 
Bouret, J.-C, Lanz, T., & Hfflier, D. J. 2005, A&A, 438, 301 
Bouret, J.-C, Lanz, T., Hfflier, D. J., et al. 2003, ApJ, 595, 1182 
Bromm, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 
Caputo, F, Marconi, M., & Ripepi, V. 1999, ApJ, 525, 784 
Cole, A. A. 1998, ApJ, 500, L137+ 
Crowther, P. A. 2000, A&A, 356, 191 
Crowther, P. A. & Dessart, L. 1998, MNRAS, 296, 622 

Crowther, P. A., Dessart, L., Hillier, D. J., Abbott, J. B., & Fullerton, A. W. 

2002a, A&A, 392, 653 
Crowther, P. A., Hillier, D. J., Evans, C. J., et al. 2002b, ApJ, 579, 774 
Crowther, P. A., Morris, P. W., & Smith, J. D. 2006, ApJ, 636, 1033 
di Benedetto, G. P. 1997, ApJ, 486, 60 

Eversberg, T., Lepine, S., & Moffat, A. F. J. 1998, ApJ, 494, 799 

Foellmi, C, Moffat, A. F. J., & Guerrero, M. A. 2003, MNRAS, 338, 360 

Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 

Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Nature, 395, 670 

Gayley, K. G. & Owocki, S. P. 1995, ApJ, 446, 801 

Genzel, R., Lutz, D., Sturm, E., et al. 1998, ApJ, 498, 579 

Grafener, G. & Hamann, W.-R. 2005, A&A, 432, 633 

Grafener, G. & Hamann, W.-R. 2008, A&A, 482, 945 

Grevesse, N. & Sauval, A. J. 1998, Space Science Reviews, 85, 161 

Groenewegen, M. A. T. 2000, A&A, 363, 901 

Hamann, W.-R. & Koesterke, L. 1998, A&A, 335, 1003 

Harries, T. J., Hilditch, R. W., & Howarth, I. D. 2003, MNRAS, 339, 157 

Herald, J. E., Hillier, D. J., & Schulte-Ladbeck, R. E. 2001, ApJ, 548, 932 

Herrero, A., Kudritzki, R. P., Vilchez, J. M., et al. 1992, A&A, 261, 209 

Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 

Hillier, D. J. 1991, A&A, 247, 455 

Hillier, D. J., Davidson, K., Ishibashi, K., & Gull, T. 2001, ApJ, 553, 837 

Hillier, D. J., Lanz, X, Heap, S. R., et al. 2003, ApJ, 588, 1039 

Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 

Hillier, D. J. & Miller, D. L. 1999, ApJ, 519, 354 

Hirschi, R., Meynet, G., & Maeder, A. 2005, A&A, 443, 581 

Hjorth, J., Sollerman, J., M0ller, P., et al. 2003, Nature, 423, 847 

Howarth, I. D. 1983, MNRAS, 203, 301 

Hunter, I., Lennon, D. J., Dufton, P. L., et al. 2008, A&A, 479, 541 

Keller, S. C. & Wood, P. R. 2006, ApJ, 642, 834 

Langer, N. 1992, A&A, 265, L17 

Lanz, T. & Hubeny, I. 2003, ApJS, 146, 417 

Lepine, S. & Moffat, A. F. J. 1999, ApJ, 514, 909 

Lepine, S. & Moffat, A. F. J. 2008, ApJ, 514, 909 

Maeder, A. 1987, A&A, 178, 159 

Marchenko, S. V., Foellmi, C, Moffat, A. F. J., et al. 2007, ApJ, 656, L77 
Marchenko, S. V., Moffat, A. F. J., Crowther, P. A., et al. 2004, MNRAS, 353, 
153 

Martayan, C, Fremat, Y., Hubert, A.-M., et al. 2007, A&A, 462, 683 
Martins, F, Hillier, D. J., Paumard, T., et al. 2008, A&A, 478, 219 
Martins, F. & Plez, B. 2006, A&A, 457, 637 

Massa, D., Fullerton, A. W., Sonneborn, G., & Hutchings, J. B. 2003, ApJ, 586, 
996 

Massey, P. & Duffy, A. S. 2001, ApJ, 550, 713 
Meynet, G. & Maeder, A. 2000, A&A, 361, 101 
Meynet, G. & Maeder, A. 2003, A&A, 404, 975 
Meynet, G. & Maeder, A. 2005, A&A, 429, 581 
Meynet, G. & Maeder, A. 2007, A&A, 464, LI 1 

Meynet, G, Maeder, A., Schaller, G., Schaerer, D., & Charbonnel, C. 1994, 
A&AS, 103, 97 

Moffat, A. F. J., Drissen, L., Lamontagne, R„ & Robert, C. 1988, ApJ, 334, 1038 
Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2007a, A&A, 465, 1003 
Mokiem, M. R., de Koter, A., Vink, J. S., et al. 2007b, A&A, 473, 603 
Najarro, F, Figer, D. F, Hillier, D. J., & Kudritzki, R. P. 2004, ApJ, 611, L105 
Nakamura, F. & Umemura, M. 1999, ApJ, 515, 239 
Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2004, A&A, 422, 675 
Oskinova, L. M., Hamann, W.-R., & Feldmeier, A. 2007, A&A, 476, 1331 
Owocki, S. P., Castor, J. I., & Rybicki, G. B. 1988, ApJ, 335, 914 
Prochaska, J. X., Bloom, J. S., Chen, H.-W., et al. 2004, ApJ, 611, 200 
Puis, J., Markova, N, Scuderi, S., et al. 2006, A&A, 454, 625 
Puis, J., Springmann, U., & Lennon, M. 2000, A&AS, 141, 23 



F. Martins et al.: WNh stars in the SMC 



13 



Runacres, M. C. & Owocki, S. P. 2002, A&A, 381, 1015 
Runacres, M. C. & Owocki, S. P. 2005, A&A, 429, 323 
Schmutz, W. 1997, A&A, 321, 268 

Schnurr, O., Moffat, A. F. J., St-Louis, N., Casoli, J., & Chene, A. . 2008, ArXiv 
e-prints, 806 

Schwarzschild, M. 1958, Structure and evolution of the stars. (Princeton, 

Princeton University Press, 1958.) 
Seaton, M. J. 1979, MNRAS, 187, 73P 

St. -Louis, N., Moffat, A. F. J., Lapointe, L., et al. 1993, ApJ, 410, 342 
St-Louis, N., Smith, L. J., Stevens, I. R., et al. 1989, A&A, 226, 249 
Steidel, C. C, Giavalisco, M., Dickinson, M., & Adelberger, K. L. 1996, AJ, 
112,352 

Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769 
Vink, J. S. & de Koter, A. 2005, A&A, 442, 587 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 

Walborn, N. R., Morrell, N. I., Ftowarth, I. D., et al. 2004, ApJ, 608, 1028 

Wellstein, S. & Langer, N. 1999, A&A, 350, 148 

Woosley, S. E. 1993, ApJ, 405, 273 

Woosley, S. E. & Bloom, J. S. 2006, ARA&A, 44, 507 

Yoon, S.-C. & Langer, N. 2005, A&A, 443, 643 

Yoon, S.-C, Langer, N., & Norman, C. 2006, A&A, 460, 199 



Appendix A: Best fits 

In this Section, we gather the figures showing the comparison 
of our best fit models with the observed spectra of our program 
stars. 
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Fig. A.l. Best fit (red solid line) of the observed spectrum (black solid line) of star SMC-WR4. The main stellar lines are indicated. 
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Fig. A.2. Best fit (red solid line) of the observed spectrum (black solid line) of star SMC-WR2. The main stellar lines are indicated. 
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Fig. A.3. Best fit (red solid line) of the observed spectrum (black solid line) of star SMC-WR1. The main stellar lines are indicated. 



